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Abstract 

A well-known problem in Malliavin calculus concerns the relation between 
the determinant of the Malliavin matrix of a random vector and the determi- 
nant of its covariance matrix. We give an explicit relation between these two 
determinants for couples of random vectors of multiple integrals. In particu- 
lar, if the multiple integrals are of the same order and this order is at most 4, 
we prove that two random variables in the same Wiener chaos either admit a 
joint density, either are proportional and that the result is not true for random 
variables in Wiener chaoses of different orders. 
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1 Introduction 



The original motivation of the Malhavin calculus was to study the existence and 
the regularity of the densities of random variables. In this research direction, the 
determinant of the so-calUed Malliavin matrix plays a crucial role. 

We give here an explicit formula that connects the determinant of the Malliavin 
matrix and the determinant of the covariance matrix of a couple of multiple stochastic 
integrals. This is related to two open problems stated in [T]. In this reference, the 
authors showed that, if F = [Fi, .., Fd) is a random vector whose components belong 
to a finite sum of Wiener chaoses, then the law of F is not absolutely continuous 
with respect to the Lebesque measure if and only if E'det A = 0. Here A denotes the 
Malliavin matrix of the vector F. In particular, they proved that a couple of multiple 
integrals of order 2 either admits a density or its components are proportional. 

They stated two open questions (Questions 6.1 and 6.2 in arXiv version): 
if C is the covariance matrix and A the Malliavin matrix of a vector of multiple 
stochastic integrals, 

• is there true that E det A > cdet C, with c > an universal constant? 

• is there true that the law of a vector of multiple integrals with components in the 

same Wiener chaos is either absolutely continuous with respect to the Lesque 
measure or its components are proportional? 

We make a first step in order to answer to these two open problems. Actually, 
we find an explicit relation that connects the two determinants. In particular, if the 
multiple integrals are of the same order and this order is at most 4, we prove that 
two random variables in the same Wiener chaos either admit a joint density, either 
are proportional. The basic idea is to write the Malliavin matrix as a sum of squares 
and to compute the dominant term of its determinant. 

We organized our paper as follows. Section 2 contains some preliminaries 
on analysis on Wiener chaos. Section 3 is devoted to express the Malliavin matrix 
as the sum of the squares of some random variables and in Section 4 we derive an 
explicit formula for the determinant of A which also involves the determinant of the 
covariance matrix. In Section 5 we discuss the existence of the joint density of a 
vector of multiple integrals. 

2 Preliminaries 

We briefly describe the tools from the analysis on Wiener space that we will need 
in our work. For complete presentations, we refer to |1] or [2]. Let H he a real and 
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separable Hilbert space and consider (W{h),h E H) an isonormal process. That is, 
{W{h),h e H) is a family of centered Gaussian random variables on the probability 
space P) such that EW{h)W{g) = {f,g)H for every h,g E H. Assume that 

the a- algebra J-" is generated by W. 

Denote, forn > 0, by Tin the nth Wiener chaos generated by W. That is. Tin is 
the vector subspace of L^(fi) generated by {Hn{W{h)), h E H, \\h\\ = 1) where Hn the 
Hermite polynomial of degree n. For any n > 1, the mapping In{h'^^) = Hn{W{h)) 
can be extended to an isometry between the Hilbert space if®" endowed with the 
norm v^t^! | | • 11//®^! and the nth Wiener chaos Tin- The random variable /«(/) is called 
the multiple Wiener Ito integral of / with respect to W. 

Consider {ej)j>i a complete orthonormal system in H and let / E H^"-, g E 
jj®m symmetric functions with n,m > 1. Then 

/= Yl ^iiv-jneji ® ••• ® ej„ (1) 
ii,..,i„>i 

and 

9= Yl /^ki,..,kA ® •• ® efc^ (2) 

where the coefficients Xi and /3j satisfy Xj^^,y...j^^„^ = Xn,..,j^ and Pk^^,^,...,k^^^^ = Pk,,..,km 
for every permutation a of the set {1, n} and for every permutation tt of the set 
{1, .., m}. Actually Xj,^„j^ = (/, ej^(g)...(g) ejj and = {g, e/,, O .. O in ([I]) 

and ([2]). Note that, throughout the paper we will use the notation (-, ■) to indicate 
the scalar product in H^'', independently of k. 

If / G if®", g E H®^ are symmetric given by ([1]), ([2]) respectively, then the 
contraction of order r of F and g is given by 

>1 

X (e^, ® .. ® ej„_,.) ® (cfc, efc„_J (3) 

for every r = 0, .., m A n. In particular f ®o g = f ® g. Note that f ®r 9 belongs to 
j^(^(m+n-2r) every r = 0, .., mAn and it is not in general symmetric. We will denote 
by f®rg the symmetrization of / ®r g- In the particular case when H = L'^{T, B, /i) 
where is a sigma-finite measure without atoms, (|3]) becomes 

dfl{ui)..dfi{Ur)f{Ui, ..,Ur,ti, . . , tn~r) 9{Ul, . . , , ^n-r+l , ..,tm+n-2r) (4) 
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An important role will be played by the following product formula for multiple 
Wiener-Ito integrals: if / G if®", g G if*^™ are symmetric, then 

mAn 

W)IM = rlC:^C:im+n-2r {f^rg) • (5) 

T- = 

We will need the concept of Malliavin derivative D with respect to W, but we will 
use only its action on Wiener chaos. In order to avoid too many details, we will just 
say that, if / is given by (II]) and In{f) denotes its multiple integral of order n with 
respect to W, then 

DIn{f) = n Xj^^„jjn-i{ej^^ ..^ejjej,. 

ji,..,j„>i 

If F, G are two random variables which are differentiable in the Malliavin sense, we 
will denote throughout the paper by C the covariance matrix and by A the Malliavin 
matrix of the random vector (F, G). That is, 

/ IIDFIII, {DF,DG)h\ 
'^-[{DF,DG)h \\DF\\l )■ 

3 The Malliavin matrix as a sum of squares 

In this section we will express the determinant of the Malliavin matrix of a random 
couple as a sum of squares of certain random variables. This will be useful in order 
to derive the exact formula for the determinant of the Malliavin matrix and its con- 
nection with the determinant of the covariance matrix for a given random vector of 
dimension 2. 

Let / G if®" and g G H'^"'- be given by ^ and ([2]) respectively, with n,m > 1. 
Let F = /„(/), G = Im{g) denote the multiple Wiener-Ito integrals of / and g with 
respect to W respectively. Then 

Inif)= Yl ^h,-,jJniej,(g) ...^ejj (6) 

and 

Imig) = Y Pk^,..,kjm (Cfci 8) .. ® tkj . (7) 

From and ([7]) we have 

DF = n Y ^ii,..,in4-i(ei2®.-®ejjej, 
ji,..,j„>i 
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and 

This implies 

\\DFfH ^n^^ \,j2,..,jnKk2,..,kJn-i {ej2 e^J 4-1 (e^^ - CfeJ 

i>l j2,-,jn>'i-kl,..,kn>l 

and 

'>1 l,32,--,jn>^lM,--,kn>l 

and 

»>1 j2,--,jn>'i- kl,..,km>l 

Let us make the following notation. For every i > 1, let 

^iJ^^Yl Yl \i2,-,3Jn-i{ej2® -^ejn) (8) 

and 

'S'i.ff = X! XI A,fc2,..,fc^/m-i (efc2 ® .. ® efc„) . (9) 

i>l i,fe2,.-,fcm>l 

We can write 

i>l l>l 

and 

i,l>l 

A key observation is that 




i,l>l \i>l J i,l>l 



(10) 



We obtained 
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Proposition 1 The determinant of the Malliavin matrix A of the vector {F, G) 
{In{f), Im{g)) can be expressed as 



il>l 



where Sij,Si^g are given by ^ and ^ respectively. 

Corollary 1 The determinant of the Malliavin matrix A of the vector {F, G) 
{In{f),Im{g)) can be expressed as 

detA = IY1 i{DF,e,){DG,ei) - {DF,ei){DG,e,)f 

i,l>l 

Proof: This comes from Proposition [I] and the relations 

S,j = {DF,e,), S,,g = {DG,e,) 

for every i >1. 



4 The determinant of the Malliavin matrix on Wiener 
chaos 

Fix n,m > 1 and /, g in iJ*^", H^"^ respectively defined by ([1]) and ([2]). Consider the 
random vector {F,G) = {In{f),Im{g)) and denote by A its Malliavin matrix and by 
C its covariance matrix. 

Let us compute i?det A. Denote, for every i,l > 1 

= ^ \h.-,u^h ® ■■® (11) 

J2vjm>l 

and 

^1,9 = ^ PiM,..,kA ® •• ® efc„. (12) 

fc2,..,fcm>l 

Clearly, for every i,l>l 

SiJ = In-l{Sij), Si^g = Im-l{Si,g)- (13) 

The following lemma plays a key role in our construction. 
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Lemma 1 /// G H^^ and g G if®'" are given by (CP and ^ respectively and Sij, Si^g 
by 07]) ■ ( fl^) respectively, then for every r = 0, .., n — 1 



nm '■ — ' 

j>i 



Proof: Consider first r = 0. Clearly, by (|3]) 



— (si,/ ® Si.g) ■ 



nm 

i>l 

The same argument applies for every r = 1, .., n — 1. Indeed, 

/ ®r+l 9 

Vil,.-Jn>l / \ki,..,km>l 

~ ^ ] ^ ] ^ ^ -^nv-,«r + l,ir + 2v- Al,..,ir + l,fcr + 2,--,fcm (^ir + 2 ® •• ® ^jn) ® (^fcr + 25 ^fcm) 

«lvi«r + l jr + 2,--,jn kr + 2,--,km 

and by ([3]) again 

E 






i>l 

*>1 \i2vjn>l / \k2,..,km>l 

nm ^ ^ ^ ] ^ ] ^ ^ \,i2,..,ir + l,jr + 2,--,jnPi,i2:--:ir + \,kr + 2,--,km 

i>l i2,..,ir + l jr + 2,--,jn kr + 2,--,km 

^ ^ ^ Kl,..,ir+lJr + 2,--Jnf^il,..,ir + l,kr+2,..,krn {^jr + 2 ^ ■■ ^ ^jn) ^ {^kr+2y ■■y^krn) 

il,..,ir + l jr + 2,--,jn kr + 2,--,km 

f ®r+l 9- 



We make a first step to compute E det A. 
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Lemma 2 Let f G H®^,g G be symmetric and denote by A the Malliavin 

matrix of the vector (F, G) = (/«(/), Im{g))- Then we have 

(n-l)A(m-l) 

£; det A = Tk 
where we denote, for k = 0, .., (m — 1) A (n — 1), 

Tk:=lYl (^--i)' (^--i)' {rn + n-2- 2k)\\\s,j^kSi,g - SlJ^kS^J' (14) 

i,l>l 

and Si J, Si^g are given by ( flil) . ( fi^) for i >1. 
Proof: By Proposition [1] and relation f|T3l) 

2detA = ^ (/„_i(sij)/m-i(sz,g) - /„_i(sij)/m-i(si,g))^ 

i,i>l 

((m-l)A(n-l) 
k\C^_^C'^_J^+n-2-2k {siJ^kSl,g - Slj^kSi,g) 
k=0 

where we used the the product formula ([5]). Consequently, from the isometry of 
multiple stochastic integrals, 



(n-l)A{m-l) 

V 

i,l>l k=0 
(n-l)A(m-l) 

= E T,. 

k=0 

■ 

For every n, m > 1 let us denote by 

(n-l)A(m-l) 
fe=l 

Remark 1 Obviously all the terms Tk above are positive, fork = 0, .., (n — 1) A(n — 1). 
We will need two more auxiliary lemmas. 



Lemma 3 Assume fi, fs G if®" and /2,/4 G H^"^ are symmetric functions. Then 
for every r = 0, .., {m — 1) A {n — 1) we have 

ifl ®n-r /s, /2 ®m-r fi) = {fl ®r f2, fs ®r fi)- 

Proof: The case r = is trivial, so assume r > 1. Without any loss of the generality, 
assume that H is L'^{T; /i) where is a sigma-finite measure without atoms. Then, 
by Q 

(/l ®n-r /3,/2 ®m-r fi) 



djjJ'iti, ..,tr) / ..,Sr) 

dl/"~''{Vi, .., Vm~r)f2{Vl, ••, fm-r, ^1, ••, tr)f4{Vl, ••, ^^m-r, Si, .., S^) 



n — 7' 



m — 7' 



(i/i" ''{Ui, ..,Un-r) / c^/i"" '"(^;i, 

^n — r J 'J^tn — r 

ifl ®r /2)(^il, ..,M„-r,fl, ..,fm-r)(/3 ®r /4)(m1, ..,M„-r,t^l, ••,^^m-r) 
(/l ®r /2,/3 ®r /4)- 



Lemma 4 Suppose fi, f/^ G if®", /2, /s G if®™ are symmetric functions. Then 

, , mAn 

~ ~ rn 'T? ' ^ — ^ 

m + n ! 

^ ' r=0 

Proof: This has been stated and proven in [3] in the case m = n. Exactly the same 
lines of the proofs apply for m ^ n. ■ 

We first compute the term Tq obtained for A; = in (|T4l) . 
Proposition 2 Let Tq be given by [T^ with k = 0. 

(n-l)A(m-l) 

To = ^nm\n\C:_,C:,_, [\\f ®. gf - \\f ®r+i gf] • 

r=0 
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Proof: From f irij) . 

To = ^{m + n-2)\^\\sij^si^g- Sij^SiJl"^ 

i,l>l 
i,l>l 

Let us apply Lemma H] to compute these norms and scalar products. We 
obtain, by letting /i = Sij = and /2 = si^g = /a (note that Sij, si^g are symmetric 
functions in if®", ]^®m respectively) 



{m + n-2)\{sij®si^g,Sij®si^g) = {m + n - 2)\{sij§si^g, si^g§Sij) 

(n-l)A(m-l) 

= (m-l)!(n-l)! ^ C'n-lCm-l(^i,/ ®r S;,,,, S^j ®^ 5;,^) 

r=0 
(n-l)A(m-l) 

= (m-l)!(n-l)! J] ^n-iC;;_ilk,/ ®. Jp. 

Analogously, for /i = s;j = /4 and /2 = Si^g = /s in Lemma H] we get 

{m + n- 2)!(s;j(g)Si,3, ■Sij^^i.s) 

(n-l)A(m-l) 

= in- mm - iy.C:_,C:,,_,\\sij ®, S.J'. 

r=0 

Next, with /i = j, /2 = si^g, /4 = s^j, /a = 

(m + n - 2)!(sij(g)Si,g, sij^Si^g) 
= {m + n-2)\{sij^SLg,Si^g^Sij) 

(n-l)A(m-l) 

^ (n - l)!(m - l)!C;_iC;,_i(si,/ 0^ sij si,g). 

r=0 
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Then 



+ n - 2)! ^ \\sij^si^g - sij^SiJ 



■ Z>1 



(n-l)A(m-l) 

(r^-i)!(m-i)!c;_ic;;„i 



r=0 



®r Sj^glP + \\Sij ®r SlJ\^ - 2{Sij ®r Si^g, Sij (g)^ S/,, 

i,l>l 

{n-l)A(m-l) 

r=0 
{n-l)A(m-l) 

2 ^ (n-l)!(m-l)!C;_iC;_, 



ii>l 



r=0 



X 



U,/>i 

{n~l)A(m-l) 



i>l 



/>1 



2 Yl (n-l)!(m-l)!C;_,C;_, 



r=0 



II 5^./ ®r Sl^gf - II J] S^J 0r Si J 



.(16) 



i>l 



Notice that, by Lemma [H for every r = 0, .., t?, — 1 



(17) 



i>l 



We apply now Lemma [3] and we get 

^ ||Si,/ (gr Si,g||^ = ^(SjJ (gr SZ,g,SiJ (gr 
i,/>l i,l>l 

i,«>l 

= C^i^iJ ®n~l~r SiJ, ^ Sl^g ®m-r-l 

and by Lemma [1] and Lemma [3l this equals 



i,l>l 



= nW||/®,^?f. 



(18) 
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By replacing ( !T7|) and ( !T8|) in (fT6|) we obtain 



To = + ||sij®sz,g-s«j®Si,3 



i,/>l 

(n-l)A(m-l) 



r=0 



As a consequence of the above proof, we obtain 

Corollary 2 For every r = 0, .., (m — 1) A (n — 1) anc? if Sij, Si^g are given by / fli]) . 
( fT^) . /ioWs that 

n-l 

n^m^ J2 CL^lC'm-l [11/ gf - 11/ ®r+l <?f ] = llSi,/®^,,^ - Slj^Si,gf . 

r=0 i,l>l 

As a consequence, for every r = 0, .., (m — 1) A (n — 1) we have 

n-l 

E^--i^™-iOl/®'-^?f -Il/®r+i^?|p] >0. (19) 

r=0 

Proof: It is a consequence of the proof of Proposition [21 ■ 
Let us state the main results of this section. 

Theorem 1 Let f G H®^,g G H'^"^(n,m > 1) be symmetric and denote by A the 
Malliavin matrix of the vector (F, G) = {In{f), Im{g))- Then 

(n-l)A(m-l) 

detA= mnm\nlC'^_^C'^_^[\\f^r9\\^-\\f(S)r+ig\\^]+Rn,m 

r=0 

where for every n,m > 1, Rn,m is given by l[T5\) . Note that Rn,m > for every 
n,m > 1. 

Proof: It follows from Proposition [2] and Lemma [2j ■ 

In the case when the two multiple integrals live in the same Wiener chaos, we 
have a nicer expression. 
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Theorem 2 Under the same assumptions as in TheoremU\ but with m = n, we have 
detA = m'detC+{mm\)' {{C^-if - iC'^'-i?) {\\f - \\f ®n-r 9f)+Rrn 

r=l 

with Rm given by 173]) . Here [x] denotes the integer part of x. 

Proof: Suppose n < m and that m is odd. The case m even is similar. From 
Theorem [1] we have 



det A = {mm\ 



\2 



(m— 1) (m— 1) 

E \\f ®r gf - Yl {ci^-iT\\f^r^i9\\' 

r=0 r=0 

"m — 1 , ^ 

(m~l) 

r=0 T — "1-1 

' ~ 2 

m— 1 

m— 1 2 



"i-i r=0 
2 



m— 1 
2 



r=0 

m — 1 

+ (mm!)2 E (^--i)' 01/ ®n-. - 11/ ®. ] 



r=l 



where we made the change of index r' = n — 1 — r in the second and third sum above. 
FinaUy, noticing that for r = we have 

m'm\' (CrO' 01/ ®o gf - 11/ ®n = det C 
we obtain the conclusion. ■ 

Example 1 Suppose m = n = 2. Then 

detA = 16[\\f^gf-\\f^2gf]+R2 
= 4detC + i?2. 

We retrieve the formula in with R2 = 32 (||/ 0i gW^ — ||/®ifi'||^) . 
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Assume m = n = 3. Then 

det A = 9 X 36 [(||/ gf - ||/ 8)3 gf) + 9 x 36 x {{Clf - l) (||/ 8)1 gf - ||/ ®2 (^f )] + R3 
= 9 det C + 9 X 36 X 3 (||/ ®i gf - ||/ ®2 gf) + Rs- 

Suppose m = n = 4. Then 

det A = 16detC + 16 X 4! X 4! ((Cg^)^ - 1) (ll/Oi c/f - ||/®3^f ) +i?4. 



5 Densities of vectors of multiple integrals 

Let us discuss when a couple of multiple stochastic integrals has a law which is 
absolutely continuous with respect to the Lebesque measure. The situations when 
the components of the vector are in the same chaos of in chaoses of different orders 
need to be separated. 

Let us first discuss the case of variables in the same chaos. In order better 
understand the relation between det A and det C we need more information on the 
terms Rm in Theorem [2l It is actually possible to compute the last term Tm-i in 

m. 



Proposition 3 Suppose m = n and let Tm-i he the term obtained in ( fl^ for r 
m — 1. Then 

T^_i = m^m\^ [\\f (S)m-igf - {f ^1 g.g®i /)] • 



Proof: From (iMll . 

i,l>l 

= \ Si^g - Sij S 

i,l>l 



i,l>l 



(m- 1)!^ 



.i.l>l 



i,l>l 



{SiJ ® SiJ, Sl^g (g) Sl^g) - Y {SiJ ® Si^g, S^g (g) Slj) 

J,1>1 i,l>l 

(m - 1)! W [(/ 01 /, g(S)ig)- (/ ®i g, g ®i /)] 
m^m!^ [11/ g\\^ - (/ ®i g, g ®i /)] 
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where we applied Lemmas [3] and [TJ ■. 

We first answer the open problem 6.2 in [T] for chaoses of order lesser than 

five. 



Theorem 3 Let m < 4 and let f.,g& jj<sm symmetric. Then the random vector 
{F,G) = {Im{f),Im{g)) does not admit a density if and only if 

det C = 0. 

In other words, the vector [F, G) does not admit a density if and only if its components 
are proportional. 

Proof: The case m = n = 1 is obvious and the case m = n = 2 follows from fT] (it 
also follows from Example [1]). Suppose m = n = 3. Then 

det A = 9detC + 9 X 36 X ((C2^)2- 1) [||/®i^f- ||/®2^||2] 

+9 X 36 [11/ ®2 - if ®2 g, g ®2 /)] + R's 

where R'^ is the term with A; = 1 in f |T4|) . Using (/ 01 g, g f) = (/ ®2 g, g ®2 /) 
(Lemma [3]) we get 

det A = 9detC + 9 X 36 X 3 [11/01 c/f - c/,^®i /)] 

-9 X 36 X 2 [11/ ®2 gf - (/ 8)2 g, g ®2 /)] + R's- 

Suppose det A = 0. Then Tq,Ti,T2 from ([H]) vanish. In particular T2 = in (HM 
and so 

||/®2^ir- (/ ®2^,^®2 /) =0. 

This implies 

9 det C + 9 X 36 X 3 [11/ ®i gf - (/ ®i g, g ®i /)] = 

and therefore det C = because \\f ^1 gW^ — {f ®i g, g ®i f) is positive by Cauchy- 
Schwarz. 

Suppose m = n = 4. 

detA = 16detC + 16x4P((Ci)2-l) [||/®i(7f -||/®3^f] 

+16 X A? [11/ ®3 g\? - if ®3 g, g ®3 /)] + K 
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where i?4 is the sum of terms obtained for k = 1 and = 2 in f|T^ . Since (/®3 g, g®3 
f) = if ®i 9, 9 ®i /) (Lemma [3]) we get 

detA = 16detC + 16xA\mCl)'-l)[\\f^^gf-{f®^g,g0,f)] 
-16 X 4!^ ((^3^)2 - 2) [11/ ®3 - (/ ®3 ^7, ^ ®3 /)] + 

Assume detA = 0. Then in particular T3 from [14]) vanishes. So 

\\f^3 9f-{f ®3t/,^®3 /) =0 

and this imphes det C = 0. 

Remark 2 For m = n > 5, we have 

detA = 25detC + 25 X 5!^ ((Cl)2 - 1) [||/®i^f - ||/®4^f] 
+25 X 5!^ ((QY - 1) [|l/®2<7f -||/®3^7f] 
+25 X 5!^ [11/ 04 - (/ ®4 ^, ^7 ®4 /)] + i?5 

//det A = then, since T^^ vanishes, we get that ||/ (8)4 g'P — (/ ®4 g-,g®i /) vanishes. 
But this is not enough. We need some additional information in order to handle the 
difference ||/ ®2 g\\^ " 11/ ®3 gW^- One possibility is to look to the terms Ti,T2,T3 
in l[T4\ ) but these terms cannot be written in a closed form, since they involve more 
complicated contractions (some "contractions of contractions"). 

Let us finish by some comments concerning the case of variables in chaoses 
of different orders. Consider {F,G) = {In{f),Im{g)) with n ^ m. First, let us note 
that i?det A = does not imply det C = 0. This can be viewed by considering the 
following example. 

Example 2 Take F = hif) and G = hih^^) where \\h\\ = 1. In this case 

detC = 2 and detA = 0. 
One can also choose F = In{h'^"') and G = Im{h®^) with n and ||/i|| = 1. 

In the case {In{f), h{g)) there is only one term in f|T4|) obtained for k = 0. It 

reads 

To = nn\ [\\f0,gf-\\f^,gf]. 

and therefore the condition for the existence of the joint density is ||/ (8)2 fi'P — 11/ ®i 

gr>o. 

The case (/„(/), 12(5')) is more complicated and needs new ideas in order to 
obtain the if and only if condition for the existence of the density of the vector. Even 
the "last term" in (fT4|) (that is, the term obtained for /c = (m — 1) A (n — 1) cannot 
be written is a nice form. 
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